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We have studied the diffuse scattering in the reiaxor Pb(Mgi/3Nb2/3)03 (PMN) using tripie-axis 
neutron scattering techniques. The diffuse scattering first appears around the Burns temperature 
Td ~ 620 K, indicating that its origin iies within the poiar nanoregions (PNR). Whiie the reiative 
intensities of the diffuse scattering around (101), (200), and (300) are consistent with those previousiy 
reported by Valchrushev et al, they are, surprisingly, entirely different from those of the lowest- 
energy transverse optic (TO) phonon. This observation led Naberezhnov et al. to claim that this 
TO mode could not be the ferroelectric soft mode. However, a recent neutron study by Gehring 
et al. has unambiguously shown that the lowest-energy TO mode does soften on coofing, and that 
the relative intensities are similar to those of PbTiOs. If the diffuse scattering in PMN originates 
from the condensation of a soft TO mode, then the atomic displacements of the PNR must satisfy 
the center of mass condition. But, the atomic displacements determined from diffuse scattering 
intensities do not fulfill this condition. To resolve this contradiction, we propose a simple model in 
which the total atomic displacement consists of two components, 5cm and 5 shift- Scm is created by 
the soft mode condensation, and thus satisfies the center of mass condition. On the other hand, 5 shift 
represents a uniform displacement of the PNR along their polar direction relative to the surrounding 
(unpolarized) cubic matrix. Within the framework of this new model, we can successfully describe 
the neutron diffuse scattering intensities observed in PMN. 



I. INTRODUCTION 

Relaxors exhibit a broad maximum in the tempera- 
ture dependence of the dielectric constant and a signif- 
icant frequency dependence of the temperature Tmax at 
which this maximum occurs. Although these peculiar- 
ities are often referred to as a "diffuse" phase transi- 
tion, no macroscopic phase transition into a ferroelec- 
tric state occurs at Tmax- Of particular interest are the 
lead-oxide class of relaxors. This class displays extraor- 
dinary piezoelectric properties which persist over a wide 
temperature range, thus presenting great appeal for de- 
vice applications. Most of the lead-oxide relaxors are 
classified as B-site complex perovskite compounds, in 
which an average B-site valence of 4-1- is realized by a 
random occupancy of two different valence cations in 
a fixed ratio. The prototypical lead-oxide relaxor sys- 
tem is Pb2+(Mg^+ Nb^^3)03~ (PMN) which, along with 

Pb(Zni/3Nb2/3)03 (PZN), shows an enormous increase 
in piezoelectric character when doped with PbTiOa. In 
spite of a decade of research, however, arguments over 
what intrinsic mechanism drives the diffuse transition re- 
main unsettled. 

One of the most important concepts related to the mi- 
croscopic properties of relaxors is that of the so-called 
"polar nanoregions" (PNR), the first experimental eVj^ 
idence for which was obtained by Burns and Dacol.El 
Through measurements of the optic index of refrac- 
tion n{T) on single crystal specimens of several dis- 
ordered ferroelectric and relaxor compounds, including 
both PMN and PZN, they observed that n{T) devi- 



ates from a linear temperature dependence at a tem- 
perature Td (600 - 650 K for PMN) far above T^ax 
(sa 265 K at 1 kHz for PMN). They proposed a model in 
which this unexpected high-temperature deviation arises 
from small, randomly oriented, very local regions of non- 
reversible polarization (the PNR) that begin to appear 
within the otherwise non-polar crystal structure below 
Td, which is often called the Burns temperature. Re- 
cent neutron inelastic scattering studies on PZN and PZN 
doped with 8% PbTiOs (PZN-8%PT) in their respective 
cubic phases at 500 K have shown that the lowest-energy 
transverse optic (TO) phonon modes are overdamped for 
reduced wave vectors q less than a characteristic wave 
vector qjj]f 0.2 A~^, but underdamped otherwiseH It 
is now believed that this damping is caused by the PNR 
because they couple strongly to the polar nature of the 
TO modes. ^ 

Vakhrushev et aln carried out neutron diffraction stud- 
ies on PMN and observed strong diffuse scattering which 
remains even at 500 K. They found that the g- width of 
the diffuse peak, which is inversely proportional to the 
correlation length, is temperature dependent, and that 
both the Bragg and diffuse peaks exhibit history depen- 
dent effects in the field-cooled (FC) and zero field-cooled 
(ZFC) regimes, similar to that of a typical spin glass. The 
transition to this glass-like state occurs around 230 K, 
which is slightly below Tmax- The diffuse scattering in 
PMN was subsequently measured around 16 reciprocal 
lattice points, from which the directions and relativia 
magnitudes of the ionic displacements were determined.tl 
The observed diffuse scattering is broader along the di- 
rection transverse to the scattering vector Q (Q±) than 
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FIG. 1: (a) Temperatum dependence of the unit cell volume, 
measured by Zhao et alB (b) Temperature dependence of the 
neutron diffuse scattering intensity (central peak) at Q = 
(2,2,0.95), measured by Naberezhnov et alU 



it is along the longitudinal direction (Q||), a feature that 
is consistent with scattering from ferroelectric fluctua- 
tions. Bonneau et alB measured both x-ray and neutron 
powder diffraction from PMN and reported that addi- 
tional diffuse scattering appears in the tails of the Bragg 
peaks below about 600 K, Neutron powder diffraction 
data taken by Zhao et alQ show a marked break in the 
linear temperature dependence of the cubic unit cell vol- 
ume, also around 600 K, andjs shown in Fig. |l|(a). More 
recently, Naberezhnov et aln have confirmed that the 
diffuse scattering indeed becomes visible below 650 K, as 
shown in Fig. |^(b). These results strongly imply that the 
diffuse scattering in PMN results from the formation of 
the PNR at T^. 

Naberezhnov et aM have also studied the lattice dy- 
namics of PMN using neutron inelastic scattering meth- 
ods between 300 and 900 K. They calculated dynami- 
cal structure factors using the room temperature diffuse 
scattering results of Vakhrushcv et alB, and concluded 
that (221) was the best zone in which to look for the 
ferroelectric soft mode phonon. However, they observed 
no underdamped soft TO phonon in the vicinity of (221) 
even at 900 K. Instead, they observed well-shaped peaks 
from the transverse acoustic (TA) and the lowest-lying 
TO modes near (220), although nearly no diffuse scat- 



tering was observed around (220). They therefore con- 
cluded that the observed lowest-lying TO mode could not 
be identified with the ferroelectric soft mode because it 
exhibited a structure factor that was absolutely incon- 
sistent with that expected from the ferroelectric diffuse 
scattering peak intensities. However, a recent neutron 
inelastic scattering study by Gehring et alE has unam- 
biguously shown that the lowest-energy zone center TO 
mode does soften significantly on cooling from 1100 K to 
Td, below which it becomes overdamped, and that the 
relative phonon intensities measured in different zones 
are very similar to those found in PbTiOs, a prototypi- 
cal (displacive) ferroelectric system. 

It is thus clear that an important discrepancy exists 
between the relative intensities of the soft TO mode and 
those of the diffuse scattering in PMN. Since each is re- 
lated to either dynamic or static atomic displacements 
toward the same ferroelectric state, the relative intensi- 
ties should be consistent. The aim of this paper is to 
determine how the diffuse scattering in PMN connects to 
the soft TO mode, and thereby resolve this discrepancy. 
To study the diffuse scattering in PMN in detail, we have 
carried out neutron diffraction measurements from 200 
to 700 K. We have concentrated on the diffuse scatter- 
ing around the three Bragg reflections (101), (200), and 
(300), and measured the g-profiles and the temperature 
dependence. These results, as well as the experimental 
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conditions, are summarized in Section In Section 
we describe our phonon structure factor calculations us- 
ing the atomic displacements determined from the diffuse 
scattering intensities, which helps to clarify the incon- 
sistency between the relative intensities of the soft TO 
mode and of the diffuse scattering. We believe this dis- 
crepancy can only be resolved if one assumes the simple 
model in which the PNR are shifted along their polar di- 
rection relative to the surrounding cubic matrix. Finally, 
we compare our experimental results with those from pre- 
vious neutron and x-ray scattering studies, and discuss 
one possible origin of this new concept, which we call a 
"phase-shifted condensed soft mode," in Section IV. 



II. DIFFUSE SCATTERING 

The neutron diffuse scattering data presented here 
were obtained on the BT9 triple-axis spectrometer lo- 
cated at the NIST Center for Neutron Research. The 
data were taken with fixed incident neutron energy Ei = 
14.7 meV (A = 2.36 A), and with horizontal beam coUi- 
mations 40'-46'-S-40'-80' or 40'-10'-S-10'-80' ("S" = sam- 
ple) . Single crystals of PMN were grown by a top-seeded 
solution growth technique using PbO as flux. The growth 
conditions were determined based on the pseuda-binary 
phase diagram established from PMN and PbO. El An as- 
grown single crystal with a nearly half-cubic morphol- 
ogy, having a volume of 0.4 cm^ and a weight of 3.25 g, 
was used for the diffuse scattering measurements. The 
crystal exhibits three naturally grown {10 0} cubic facets, 
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FIG. 2: (a) Diffuse scattering profiles around (101) along [101] 
at 300, 400, 500, 600, and 690 K. Curves are fits to a com- 
bination of Gaussian (Bragg peak) and Lorentzian (diffuse 
scattering) lineshapes. (b) Temperature dependence of the 
diffuse scattering intensities at (1.20, 0, 0.80), (1.15, 0, 0.85) 
and (1.10, 0, 0.90). The diffuse scattering clearly begins in- 
creasing around Td- 



and was mounted on a boron nitride post using tanta- 
lum wire with one of the facets facing verticaUy, and then 
attached to the cold-head of a high-temperature closed- 
cycle helium refrigerator. This orientation gave access to 
reflections of the form (hOl) in the scattering plane. The 
lattice constant of PMN is a = 4.04 A at room temper- 
ature, thus 1 rlu (reciprocal lattice unit) corresponds to 
27r/a = 1.553 k-\ 

Figure |^(a) shows peak profiles of the transverse diffuse 
scattering in the vicinity of (101). Note that the y-axis 
is displayed on a logarithmic scale. At 690 K, there is no 
diffuse scattering. The peak profile is well-described by 
a Gaussian function with a full-width at half-maximum 
(FWHM) corresponding to that of the instrumental q- 
resolution at the (101) Bragg peak. With decreasing tem- 
perature, a weak signal starts to emerge out of the (101) 
Bragg peak, growing more rapidly as the temperature is 
lowered. The observed diffuse scattering can be nicely fit 
with a Lorentzian function as shown in Fig. ^(a). The 
temperature dependence of the diffuse scattering is more 
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FIG. 3: Profiles of the diffuse scattering peaks at (a) (1 — 
q,0,l + q), (b) (2,0,g), and (c) (3, 0, g) at T = 300 and 370 K. 



clearly presented in Fig. @(b). The diffuse scattering ap- 
pears below about 600 — 650 K, consistent with the Burns 
temperature Td ^ 620 K for PMN, and increases almost 
exponentially with decreasing temperature. 

The diffuse scattering peak profiles were more closely 
examined by employing a tighter coUimation of 40'-10'-S- 
10'-80'. Figure H depicts the transverse diffuse scattering 
profiles at (101), (200), and (300) at 300 and 370 K. 
Although the difference in the diffuse scattering inten- 
sity between these two temperatures is not very large, 
as is apparent from Fig. 0, it is nevertheless discernable 
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FIG. 4: Intensity contours of the diffuse scattering around 
(101) at (a) 270 K and (b) 370 K. (c) Difference between (a) 
and (b). 



at (101). A significantly larger difference is observed 
at (300), which indicates that the diffuse scattering is 
stronger at (300) than at (101). No difference in inten- 
sity at finite q near (200) was found, which suggests the 
diffuse scattering is very weak at (200). We speculate 
that the tail-like feature near g = for (200) is due to 
contamination from the TA phonon, which is quite strong 
near (200), as some TA mode scattering will inevitably 
spill into the elastic profile because of the imperfect in- 
strumental energy resolution. The relative intensities of 
the (101), (200), and (300) diffuse scattering peaks are 
consistent with those reported in the previous neutron 
diffraction measurement by Vakhrushev et ala We have 
also carried out a detailed survey of the (101) diffuse 
scattering, the results of which are shown in Fig. ^. The 
diffuse scattering is highly elongated along the direction 
transverse to the scattering vector, which is also consis- 
tent with the previous report Similar neutron diffuse 
scattering results at (101) and (200) were recently ob- 
served by Koo et al. in ;i-aingle crystal sample of PMN 
doped with 20% PbTi03.y 
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FIG. 5: Energy spectra of the transverse acoustic and op- 
tic phonons at 690 K, i.e., above the Burns temperature T^, 
around (a) (200) and (b) (300). 



III. FERROELECTRIC PHONON 

Gehring et aZ.i have recently reported that the lowest- 
energy zone center TO mode in PMN softens on cool- 
ing from 1100 K to Td, and that the relative intensi- 
ties measured in different zones are very similar to those 
found in PbTiOa. This result sharply contradicts the 
fact that the relative intensities of the diffuse scattering 
in different zones are entirely different from those of the 
lowest-energy transverse optic (TO) jphonons, an obser- 
vation which led Naberezhnov et alu to claim that the 
lowest-energy TO mode could not be the ferroelectric soft 
mode. To understand how PMN undergoes the diffuse 
transition through the formation of the PNR, this con- 
tradiction must be resolved. In this section, we revisit 
the soft mode of PMN and introduce a new concept, the 
"phase-shifted condensed soft mode," as a microscopic 
description of PNR formation. 

As shown in the previous section, we observe no diffuse 
scattering around (200), consistent with previous reports. 
However, TO phonons are clearly observed around both 
(200) and (300) as shown in Fig. |, which shows two 
constant-Q scans at 690 K (i.e., above the Burns tem- 
perature Td) at (2,0,-0.2) and (3,0,-0.2). The solid 
curves are fits to two Lorentzian functions of i?, the neu- 
tron energy transfer, convolved with the proper instru- 
mental resolution function. From the fitting results, the 
ratio of the two TO phonon intensities is 

^^«g = 1.24 ±0.20. (1) 

|i^ofc.(300)|2 ^ ^ 

According to the pioneering work by Harada et aZ.0, 
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FIG. 6: (a) Calculated intensities of the TO phonon mode at 
(101), (200) and (300) as a function of the ratio S between the 
Si (Slater) and 5*2 (Last) modes, (b) Calculated intensities of 
the TO phonon mode at (101), (200) and (300) as a function of 
the phase shift of the PNR along their polar direction relative 
to the surrounding cubic matrix. 



which describes the determination of the normal mode 
vibrational displacements in perovskites from measured 
phonon intensities, the relative intensities of phonons 
are determined mostly by the ratio of two dominant 
modes, i.e., the Slater mode and the Last mode. In 
the Slater mode, the oxygen and Mg/Nb (MN) atoms 
vibrate in opposition while the Pb atoms remain station- 
ary. The Last mode corresponds to opposing motions of 
the (Mg/Nb)06 octahedra and the Pb atoms. In both 
modes, the three oxygen atoms in the unit cell move as 
a rigid unit. We define the ratio of the two modes as 
S = S2/S1, where Si and S2 represent the contribution 
from the Slater mode (si) and Last mode (S2), respec- 
tively. 

Let us first determine the ratio S. Figure |^(a) shows 
the inelastic structure factors |Fi„eip at (101), (200) 
and (300) as a function of S, which-are calculated using 
the formula given by Harada et alSA In this calculation, 
we assume that the atomic displacements Si and S2 are 
parallel to the [100] direction for the sake of simplicity. 
Combining Eqs.(l) and (5) of Ref. 11, we have obtained 



\Finei\ for the lowest branch at the zone center (q — 0) 
using constant- Q scans as follows: 



\F„ 



(2) 



rocal lattice vector (G) and the phonon wave number 
(q = 0), and the normalized displacement vector for the 
K-th atom, in our condition is 



(3) 



From Table 1 of Ref. |Tl|, the atomic displacements along 
the [100] direction are 



Co 



~Sk' 
-k + S 
l + S, 



(4) 



where the center of mass condition requires k = 
3M(0)/M(MN) = 0.686 and k' = [Af(MN) + 
3A//(0)]/M(Pb) = 0.596. The calculated values are nor- 
malized by the sum of \Finei\'^ over all reflections that 
can be reached with an incident neutron energy Ei — 
14.7 meV: (100), (110), (200), (210), (220), (300), (310) 
and their equivalent. Debye- Waller factors, exp {—Wk,) — 
exp(— BK|Q/47r|^), are calculated ysing Bpi, — 0.915 
and Bnm = 0.483 and Bq = 1.02113 Figure |(a) indi- 
cates that the experimentally observed intensity ratio of 
|Fofc,(200)|V|Fob^(300)|2 = 1.24 obtained in Eq. (|l|) is 
realized between = 1.0 and 1.5. 

If the diffuse scattering of PMN originates from the 
condensation of the soft TO mode, then the above for- 
mula should be also applicable to the diffuse scattering 
intensities and that the atomic displacements determined 
from the diffuse scattering intensities must also satisfy 
the center of mass condition: 



^J(k)M(k) = 0, 



(5) 



where 5{k) and M(k) are the displacement and the 
atomic mass of the K-th atom, respectively. However, 
the values S{k) determined from the diffuse scattering 
intensities by Vakhrushev et aln, 



S(Ph) = 1.00, 
(5(MN) = 0.18, 
(5(0) = -0.64, 



Af(Pb) = 207, 
Af(MN) = 70, 
M(0) = 16, 



(6) 



do not satisfy the center of mass condition. This con- 
tradiction is a natural corollary to the observed inconsis- 
tency between the relative intensities of the diffuse and 
the soft TO mode scattering. 

To resolve this discrepancy, we propose a simple model 
in which the atomic displacements consist of two compo- 
nents, S{k) — 6cm{h) + Sshift- In this model, Scm is 
induced by the condensation of the soft mode, and thus 
satisfies the center of mass condition. On the other hand, 
Sshift shifts the PNR along their polar directions rela- 
tive to the surrounding cubic matrix, which is, in other 
words, a phase shift. By applying this model, the S{k) 
values listed in Eq. (o) can be separated into two parts; 



where the scattering vector (Q) is the sum of the recip- 



ScMiPh) = 0.42, 
Scm(MN) = -0.40, 
ScMiO) = -1.22, 



(7) 



and 



^shift 



0.58. 



(8) 



Note that the 5cm values are uniquely determined from 
the center of mass condition and can be further decom- 
posed into contributions from the Slater mode and the 
Last mode. The corresponding ratio S between the Slater 
mode and Last mode contributions is 1.5, which is consis- 
tent with the value obtained in Fig. ^(a). Note that the 
phase shift 5 shift becomes effective only when the PNR 
condense out from the soft TO mode. Hence the phonon 
intensities are determined only by the 5cm values. 

Now that we have confirmed that the 5cm values de- 
rived from the diffuse scattering intensities are consistent 
with the inelastic scattering intensities of the soft TO 
mode, we can calculate how the diffuse scattering intensi- 
ties depend on the phase factor 5shi ft with the constraint 
S" = 1.5 using the 5cm values listed in Eq. (0). The re- 
sults are shown in Fig. ^(b). Note that the values at 
5 shift = are identical to those at S" = 1.5 in Fig. ||(a) as 
expected. At 5shift = 0.58, which was derived above to 
satisfy the center of mass condition, the diffuse intensity 
around (200) approaches zero, in agreement with exper- 
iment. Thus, our "phase-shifted condensed soft mode" 
model consistently explains both the diffuse scattering 
and the soft mode intensities from the same atomic dis- 
placements. 

Lastly, we consider the phase shift from the crystal 
lattice point of view. Figure |^ depicts the atomic dis- 
placements schematically with three different phase shift 
values. Figure |^(a) shows the displacements of Eq. (|^) 
without the phase shift, which satisfies the center of mass 
condition and gives the same relative intensities for both 
the diffuse scattering and the soft TO mode. However, 
once all the atoms undergo an additional shift along the 
same direction by the same amount 5 shift = +0.6, the 
center of mass condition is no longer satisfied, resulting 
in the accidental disappearance of the (200) diffuse neu- 
tron scattering intensity. Figure 7|(c), which depicts a 
uniform shift of 5 shift = —0.6, is shown just for compar- 
ison. 



IV. DISCUSSION 

In the preceding sections we have demonstrated how 
the inconsistency between the diffuse scattering and the 
soft TO mode intensities is explained naturally within 
the framework of a "phase-shifted condensed soft mode." 
Our data on the diffuse scattering cross-sections in three 
different Brillouin zones along with our measurements of 
the soft TO mode phonons in two different zones serve 
as a confirmation of the detailed measurements already 
published by other groups. eIeIEtlI But while our data are 
in agreement, our interpretations differ entirely. The sim- 
plicity of our model argues in its favor. Other than the 
standard picture of condensed local ferroelectric fluctua- 
tions created by a soft phonon, there is only one added 
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feature, that being a uniform phase shift of the PNR rel- 
ative to the surrounding cubic matrix. We believe the 
simpler solution is the correct one. 

The process of separating the observed atomic dis- 
placements into optic mode displacements 5cm and a 
single phase shift parameter SsM ft is straightforward and 
unique, and requires no adjustable parameters. We re- 
alized this fact a long time ago. We did not recognized, 
however, the significance of the derived center of mass 
structure. Only recently, when the true soft mode of 
PMN was observed directly by neutron scattering at high 
temperaturejj did we recognize the identity of the diffuse 
structure as a condensed soft mode. Then the idea of the 
phase shift followed naturally. But one important ques- 
tion still remains: what is the microscopic origin of this 
phase shift? It should be emphasized that the direction 
of the phase shift is not arbitrary as is clear from Fig. |^. 
Our model requires that the phase shift be parallel, not 
antiparallel, to the polarization resulting from the atomic 
displacement. Under zero external electric field, as is the 
case here, the polarization direction of a single PNR can 
be parallel to aw symmetric direction of (111), i.e., [Ill], 
[ill],- ■ •, [iil].E£l Therefore, the phase-shifted directions 
of the PNR are microscopically aligned to their polar- 
izations, although they appear macroscopically random. 
Then, why is the phase shift of the PNR parallel to their 



7 



polarization? 

We speculate that microscopic inhomogeneities in the 
site-occupancy of the Mg^+ and Nb^+ cations produce a 
local electric field gradient that determines the direction 
of polarization of the PNR. In other words, a nanoregion 
exposed to a local field gradient becomes a polar nanore- 
gion with a polarization parallel to the gradient. Once 
a PNR is formed, its polarization starts to interact with 
the local field gradient, and may result in shifting the 
PNR in order to compensate for the Coulomb interac- 
tion with the lattice distortion due to the phase shift. 
The phase shift in this case is parallel to the local field 
gradient, and thus to the polarization of the PNR. Since 
the microscopic inhomogeneities associated with the spa- 
tial distribution of the Mg^+ and Nb^+ cations causes 
random charged domains, the local field gradient is ran- 
dom. Therefore, the polarizations of the PNR should 
also be random. However, as explained above, the phase 
shift of a particular PNR is always parallel to its polar- 
ization. These descriptions are still qualitative, and need 
to be tested experimentally. It is particularly important 
to study the microscopic properties of the PNR, such as 
their correlation lengths and textures as a function of 
temperature as well as electric field. We also need to un- 
derstand the chemical inhomogeneity of Mg^+ and Nb^+ 
microscopically. 

Finally, our model should help to resolve some_Q£lli£ 
conflicting x-ray diffuse scattering interpretations £3ll3'E£l 
The diffuse intensity around (200) is only accidentally 
zero when measured with neutrons, but it should be non- 
zero when measured with x-rays because of the differ- 
ences between the neutron nuclear scattering lengths and 
the x-ray atomic scattering factors. 

In summary, we have studied the diffuse scattering of 
the relaxor Pb(Mgi/3Nb2/3)03 over a wide temperature 
range with neutron scattering techniques. We have con- 



firmed that the relative intensities are consistent with 
previous reports and that the diffuse intensity starts in- 
creasing below the Burns temperature Td- We have re- 
visited the soft transverse-optic mode, the existence of 
which has been recently confirmed, and examined the in- 
consistency between the soft phonon and the diffuse in- 
tensities. We have proposed a new concept, the "phase- 
shifted condensed soft mode," which naturally explains 
the inconsistency as well as the microscopic origin of the 
formation of the polar nanoregions. 
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